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C ' Abstract 

A specific algebraic realization of the Ginsparg- Wilson relation in the form 
75(75-^) + (75-^)75 = 20^^^+^ (75!))^^"^^ is discussed, where k stands for a non- 
negative integer and A; = corresponds to the commonly discussed Ginsparg- Wilson 
■ relation. From a view point of algebra, a characteristic property of our proposal 

^ . is that we have a closed algebraic relation for one unknown operator D, although 

this relation itself is obtained from the original proposal of Ginsparg and Wilson, 
O ' IbD + -D75 = 2aD^^aD, by choosing a as an operator containing D ( and thus 

. Dirac matrices). In this paper, it is shown that we can construct the operator D ex- 

\ plicitly for any value of k. We first show that the instanton-related index of all these 

<^ ' operators is identical. We then illustrate in detail a generalization of Neuberger's 

overlap Dirac operator to the case k = 1. On the basis of explicit construction, it is 
shown that the chiral symmetry breaking term becomes more irrelevent for larger 
pi 1' k in the sense of Wilsonian renormalization group. We thus have an infinite tower 

^ ■ of new lattice Dirac operators which are topologically proper, but a large enough 

lattice is required to accomodate a Dirac operator with a large value of k. 
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b : 1 Introduction 



Recent developments in the treatment of fermions in lattice gauge theory are based on a 
hermitian lattice Dirac operator 75Z) which satisfies the Ginsparg- Wilson relation[l] 

75D + D75 = 2aD-ir,D (1.1) 

where the lattice spacing a is utilized to make a dimensional consideration transparent, 
and 75 is a hermitian chiral Dirac matrix. An explicit example of the operator satisfying 
(1.1) and free of species doubling has been given by Neuberger[2]. The relation (1.1) led 

^To be precise, the general relation "f^D + Dj^ = 2aD^^aD, where a is a local operator, has been 
proposed in Ref. [1], although the authors in Ref.[l] analyzed "only the simplest case where the matrix a 
is proportional to the unit matrix in Dirac space". With this quaUfication in mind, we refer to (1.1) as 
the "ordinary Ginsparg- Wilson relation" in this paper. The original Ginsparg- Wilson relation is more 
general as stated above. 
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to an interesting analysis of the notion of index in lattice gauge theory[3]. This index 
theorem in turn led to a new form of chiral symmetry, and the chiral anomaly is obtained 
as a non-trivial Jacobian factor under this modified chiral transformation [4]. This chiral 
Jacobian is regarded as a lattice generalization of the continuum path integral[5]. The 
very detailed analyses of the lattice chiral Jacobian have been performed [6]- [8]. It is also 
possible to formulate the lattice index theorem in a manner [9] analogous to the continuum 
index theorem[10] [11]. An interesting chirality sum rule, which relates the number of zero 
modes to that of the heaviest states, has also been noticed[12]. 

In this paper we discuss a generalization of the relation (1.1), which is characterized 
by a non-negative integer k. It is shown that the explicit construction of an infinite tower 
of lattice Dirac operators which satisfy the index theorem is possible, but a large enough 
lattice is required to accomodate a Dirac operator with a large value of k. 

2 Generalized algebra and its representation 

We discuss a generalization of the algebra (1.1) to the formQ 

75(75!^) + (75/^)75 = 2a''^\j,Dy'-^' (2.1) 

where k stands for a non-negative integer and A; = corresponds to the ordinary Ginsparg- 
Wilson relation. When one defines 

H = j^aD (2.2) 

j,H + Hj, = 2H^^+' (2.3) 
T5H + T5H = (2.4) 



(2.1) is rewritten as 
or equivalently 
where we defined 

T,^^,-H''^^\ (2.5) 

Note that both of H and are hermitian operators. 

We now discuss a general representation of the algebraic relation (2.4) following the 
analysis in Appendix of Ref. [13].(In Ref.[13], the algebra was normalized as 75(75/?) -|- 
(75/^)75 = 0(75/?)^, but here we use the normalization (2.1) to simplify various expres- 
sions.) The relation (2.4) suggests that if 

H(f)n = aXn4>n, {4>n,4>n) = 1 (2.6) 



^This relation is obtained from the proposal in Ref. [l],75l? + I?75 = 2aD^^aD , by choosing a as 
an operator containing D itself (and thus Dirac matrices). From a view point of algebra, the original 
construction in [1] contains two unknown operators and one relation. In our construction, we have a 
closed algebraic relation for one unknown operator D, which allows a neat analyis of representation 
in this Section. This specific algebraically closed realization , which is characterized by a non-negative 
integer, has not been discussed in Ref.[l]. 
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with a real eigenvalue aA„ for the hermitian operator i?, then 

H{T^(t)n) = -a\niT5(Pn). (2.7) 

Namely, the eigenvalues A„ and — A„ are always paired if A„ 7^ and (Fs^^, r50„) ^ 0. 
We also note the relation, which is derived by sandwiching the relation (2.3) by 

(0n,750n) = (aA„)2'=+^ foT A„ ^ 0. (2.8) 

Consequently 

\{aX^r^'\ = |(0„,750n)| < ||<^n||||750n|| = 1- (2.9) 

Namely, all the possible eigenvalues are bounded by 

|A„| < -. (2.10) 

(X 

We thus evaluate the norm of T^(f)n 

(r50„,r50„) = (0„,(T5-i^''+')(75-^''^')0n) 

= (0„, (1 - i/^^+Ss - 75^^^'+^ + i^^(^'+^))0n) 

= [l-(aA„)W)] 

= [l-(aA„)2][l + (aA„)2 + ... + (aAX] (2-11) 

where we used (2.8). By remembering that all the eigenvalues are real, we find that 0^ is 
a "highest" state 

r^cPn = (2.12) 

only if 

[1 - (aA„)'] = (1 - aA„)(l + aA„) = (2.13) 

for the Euclidean positive definite inner product (0„,0n) = Sx- 0n(^)0n(^)- 

We thus conclude that the states 0„ with A„ = ±^ are not paired by the operation Fs^^ 

and ^ 

-f5D(l)n = ±-0„, 75</'n = ±(f>n (2.14) 

a 

respectively. These eigenvalues are in fact the maximum or minimum of the possible 
eigenvalues of H/a due to (2.10). 

As for the vanishing eigenvalues = 0, we find from (2.4) that -ff750n — 0, namely, 
H[{1 ± 75) /2]0„ = 0. We thus have 

75-D0„ = 0, 750n = 0n OT 750n = -0n- (2-15) 



To summarize the analyses so far, all the normalizable eigenstates 0„ of 'j^D = H/a are 
categorized into the following 3 classes: 
(i) n± ("zero modes"), 

75-D0n = 0, 750„ = ±0„, (2.16) 
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(ii) N± ("highest states"), 

^5D(f)n = ±-0n, 750n = ±4>n, respectively, (2-17) 
a 

(iii) "paired states" with < |An| < 1/a, 

-f^Dcpn = Xn(l>n, 75^(r50n) = -X„{T5(pn). (2.18) 

Note that r5(r50„) oc 0„ for < |An| < 1/a. 
We thus obtain the index relation [3] [4] 

n 

= i^n, r50n) + (0 ) + i^n, Ts^n) 

An=0 0<|An|<l/a |A„|=l/a 

A„=0 

= E (75 - /f''-'')0n) 

An=0 

A„=0 

= n+ — n_ = index (2-19) 

where n± stand for the number of normahzable zero modes with 750n = =t0n in the 
classification (i) above. We here used the fact that Fg^^ = for the "highest states" 
and that 0„ and Fs^^ are orthogonal to each other for < |An| < 1/a since they have 
eigenvalues with opposite signatures. 

On the other hand, the relation Tr75 — 0, which is expected to be valid in (finite) 
lattice theory, leads to ( by using (2.8)) 

= E(<^n'750n) 

n 

= E (^"'750n) + E (^'»'750n) 
A„=0 A„^0 

= n+ - n_ + E («A„)"=+' = 0. (2.20) 

In the last line of this relation, all the states except for the "highest states" with A„ = ±l/a 
cancel pairwise for A„ 7^ 0. We thus obtain a chirality sum rule[12] n+ — n^ + N+ — N_ — 
or 

n+ + N+=n- + N_ (2.21) 

where N± stand for the number of "highest states" with 750„ = ±0n in the classification 
(ii) above. These relations show that the chirality asymmetry at vanishing eigenvalues is 
balanced by the chirality asymmetry at the largest eigenvalues with |A„| = 1/a. It was 
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argued in Ref.[13] that N± states are the topological (instanton-related) excitations of the 
would-be species doublers. 

All the rij- and N± states are the eigenstates of D, D(f)n = and -D0„ = (l/a)0„, 
respectively. If one denotes the number of states in the classification (iii) above by 2Nq, 
the total number of states (the dimension of the representation) N is given by 

N = 2{n+ + N+ + No) (2.22) 

which is expected to be common to all the algebraic relations in (2.1) and to be a constant 
independent of background gauge field configurations. 

We note that all the states 0„ with < |A„| < 1/a, which appear pairwise with 
A„ = ±|A„|, can be normalized to satisfy the relations 

r50„ = [l-(aA„)^(^'=+^)]VV_„, 

750n = (aA„)^'=+Vn+[l-(aA„)^(^'=+^)]^/V-n. (2.23) 

Here (p-n stands for the eigenstate with an eigenvalue opposite to that of 4>n- These states 
(pn cannot be the eigenstates of 75 since |(0n,750n)| = K'^-^n)^'^'''^! < 1- 

We have thus established that the representation of all the algebraic relations (2.1) 
has a similar structure. In the next Section, we show that the index n+ — n_ is identical 
to all these algebraic relations if the operator j^D satisfies suitable conditions. 

3 Chiral Jacobian and the index relation 

The Euclidean path integral for a fermion is defined by 

j VijVip exp[ j TpDip] (3.1) 

where 

/ iljDilj = ^ilj{x)D{x,y)ilj{y) (3.2) 

x,y 

and the summation runs over all the points on the lattice. The relation (2.4) is re- written 
as 

TsFsTsi^ + i^Fs - (3.3) 
and thus the Euclidean action is invariant under the global "chiral" transformation [4] 

z 
w 

with an infinitesimal constant parameter e. Under this transformation, one obtains a 
Jacobian factor 

PV^W = JVi)Vil) (3.5) 
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with 

J = exp[— 2iTrer5] = exp[—2ie{n^ — (3-6) 

where we used the index relation (2.19). 

We now relate this index appearing in the Jacobian to the Pontryagin index of the 
gauge field in a smooth continuum limit by following the procedure in Ref.[9]. We start 
with 

Tr{T,f{^^)} = Tr{T,f{^^)} = - n_ (3.7) 

Namely, the index is not modified by any regulator f{x) with /(O) = 1 and f{x) rapidly 
going to zero for x oo, as can be confirmed by using (2.19). This means that you can 
use any suitable f{x) in the evaluation of the index by taking advantage of this property. 
We then consider a local version of the index 

H^^fi^^m^, ^) = tr{{l, - H''^')f{^^^^)}{x, x) (3.8) 

where trace stands for Dirac and Yang-Mills indices; Tr in (3.7) includes a sum over 
the lattice points x. A local version of the index is not sensitive to the precise boundary 
condition , and one may take an infinite volume limit of the lattice in the above expression. 

We now examine the continuum limit a ^ of the above local expression (3.8)0. We 
first observe that the term 

tr{H^'^'f{^^)} (3.9) 

goes to zero in this limit. The large eigenvalues of H = a'y^D are truncated at the value 
aM by the regulator f{x) which rapidly goes to zero for large x. In other words, the 

75^) " 

M2 



global index of the operator TrH'^^+^ f{^^^) ~ 0{aMf^+^. 



We thus examine the small a limit of 

tr{75/(^^)}. (3.10) 

The operator appearing in this expression is well regularized by the function f{x) , and 
we evaluate the above trace by using the plane wave basis to extract an explicit gauge 
field dependence. We consider a square lattice where the momentum is defined in the 
Brillouin zone 

TT 37r 
— < kn< — . 
2a - ^ 2a 

We assume that the operator D is free of species doubling; in other words, the operator 
D blows up rapidly (~ ^) for small a in the momentum region corresponding to species 
doublers. The contributions of doublers are eliminated by the regulator f{x) in the above 
expression, since 

Hl.n'^)) ~ (i)V(^) ^ (3.12) 



- TT7 < < S^- (3-11) 



•^This continuum limit corresponds to the so-called "naive" continuum limit in the context of lattice 
gauge theory. 
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for a — > if one chooses f{x) — e ^, for example. 

We thus examine the above trace in the momentum range of the physical species 



^<.,<iL. (3.13) 



We obtain the limiting a — > expression 



lim ir{75/ ( ^^^^ )} (x, x) 



lim limtr [' -^e-^^^75/(^^|^)e^'=^ 



= tr{75/(|^)} (3.14) 

where we first take the limit a with fixed in — L < A;^ < L, and then take the limit 
L oo. This procedure is justified if the integral is well convergent We also assumed 
that the operator D satisfies the following relation in the limit a — > 

De^^^'hix) }^ + I ^ - g ^)h{x) 

= i{^ + ig4){e''''-'h{x))=ip{e'''''h{x)) (3.15) 

for any fixed kn, (— ^ < < ^), and a sufficiently smooth function h{x). The function 
h{x) corresponds to the gauge potential in our case, which in turn means that the gauge 
potential Af^{x) is assumed to vary very little over the distances of the elementary lattice 
spacing. 

Our final expression (3.14) in the limit M — > oo reproduces the Pontryagin number in 
the continuum formulation 

hm tr{^,f{]/)'/M')} 

M—*oo 
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To be precise, we deal with an integral of the structure /_^\ dxfa{x) = J^" dxfa{x) + 
J^j^dxfa{x) + dxfa{x) where fa{x) depends on the parameter a. (A generalization to a 4- 
dimensional integral is straightforward.) We thus have to prove that both of lima^o J dxfa{^) 
and linio^o /_i dxfa{x) can be made arbitrarily small if one lets L to be large. A typical integral 
we encounter in lattice theory has a generic structure lima^o /^f /"a aa;+(i-cos2a2;) ]/(a m ) _ 

lima_o f^/f;' rfa;e-[sin^«- + (l-cos2a.)^]/(a^M^) ^ p-/2a ^ [sin^ ax+(l-cos 2ax)^]/(a^M^) ^ 

" ^ Je7r/2a " ^'J— 7r/2a 

lim n f'^''/^" J -[sin-ax + (l~cos2axf]/{a^M^) ^ rCT/2a ^^^-[sin''ax + (l-cos2axf]/{a^M^) ^ 

J— e7r/2a ^^^'J— e7r/2a 

liniL^oc J_j^dxe~^^^^^ and satisfies the above criterion, if one chooses the regulator f{x) = e^^: Here 
e is an arbitrary small fixed parameter, and the left-hand side of this relation stands for a conventional 
lattice calculation and the right-hand side stands for a continuum calculation. 
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= lim tr I e-''^^Ji^yM')e"^^ 

= ^^Jr J (0r75/{(iA:, + D.f/M' + Y]F,./M'} 

= ^rnJrM' J ^.l.fiHK + DjMf + ^[7^ Y]FjM'} (3.16) 

where the remaining trace stands for Dirac and Yang-Mills indices. We also used the 
relation 

ip^ = D^D^ + '-l[Y.Y]F,. (3.17) 

and the rescaling of the variable Mk^. 

By noting tr75 = tr75[7'^, 7^^] = 0, the above expression ( after expansion in powers of 
1/M) is written as (with e^^^'' = 1) 

]^^Jn,f{P'/M') = tr75^{f[7^7li^M4V(0/"(-^'^^') 

= ^^tre^-'^^F^^F^p (3.18) 

where we used 

r rl'^h 

J {2ny' ^ ^ ' 



with X — —k^k'^ > Q in our metric. 

When one combines (3.7) and (3.18), one reproduces the Atiyah-Singer index theorem 
(in continuum space) [10] [11]. We note that a local version of the index (anomaly) is 
valid for Abelian theory also. The global index (3.7) as well as a local version of the index 
(3.8) are both independent of the regulator f{x) provided [5] 

/(O) = 1, /(oo) = 0, f{x)x\,=o = f'{x)xl=^ = 0. (3.20) 

Wc have thus established that the lattice index in (3.7) for any algebraic relation in 
(2.1) is related to the Pontryagin index in a smooth continuum limit as 

n+-n_= j d^x^tre^'^^^F.^F^^ (3.21) 

by assuming the quite general properties of the basic operator D only: The basic relation 
(2.1) with hermitian 75/) and the continuum limit property (3.15) without species doubling 
in the limit a 0. This shows that the instanton-related topological property is identical 
for all the algebraic relations in (2.1), and the Jacobian factor (3.6) in fact contains the 
correct chiral anomaly. (We are implicitly assuming that the index (3.7) does not change 
in the process of taking a continuum limit.) Our result is naturally consistent with the 
calculation of chiral anomaly by different methods in [1] and [3]. 
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167r2 
1 

T67r2 



/ f"ix)xdx 
Jo 



(3.19) 
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4 Explicit example of the lattice Dirac operator with 

k=l 



Wc now discuss an explicit construction of the lattice Dirac operator which satisfies the 
generalized algebraic relation (2.1) with k = 1, though a generalization to an arbitrary k is 
straightforward as is described in Section 5 later. For this purpose, we first briefiy review 
the construction of the Neuberger's overlap Dirac operator for the ordinary Ginsparg- 
Wilson relation. 

We start with the conventional Wilson fermion operator D^r defined by 

Dw{x, y) = «7^C^(a;, y) + B{x, y) - -mo^a:,^, 

a 

B{x, y) = — YPKy - Sy+fia,xUl{x) - Sy,:c+f.aUi,{y)], 

Uf,{y) = exp[iagA^{y)], (4.1) 

where we added a constant mass term to D^y for later convenience. The parameter r 
stands for the Wilson parameter. Our matrix convention is that are anti-hermitian, 
{l^y — ~7^) and thus = ^^C^{n, m) is hermitian 

=(j^. (4.2) 

The operator D introduced by Ncubcrger[2], which satisfies the conventional Ginsparg- 
Wilson relation (1.1), has an explicit expression 

aD = \[l + 75^] = ^[1 + ] ] (4.3) 

where Dw = J^Hw is the Wilson operator defined above, and Hw is hermitian Hl^r = Hw- 
The physical meaning of this construction becomes more transparent if one considers 
(naive) near continuum configurations specified by a small a limit with the parameters 
r/a and mo /a kept finite. We can then approximate the operator D^r by [14] 

Dwc^ip + Mn (4.4) 

for each species doubler, where the mass parameters M„ stand for Mq — and one of 

2r mo , . 4r uiq 

(4,-1); (1,1) (4.5) 

a a a a 

for n = 1 ~ 15; we denoted ( multiplicity, chiral charge ) in the bracket for species 
doublers. Here we used the relation valid in the near continuum configurations for the 
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physical species, for example, 



Dw{k) = l^Y ^ + - cos ak^) 

^ 7%-— (4.6) 

in the momentum representation with vanishing gauge field. 

In a symbolic notation, one can then write the overlap Dirac operator as 

~ ^i[l + (z^ + M„) ^ ]|n)H, 
n=o2 

a^,D ^(-l)"75^[l + (i^ + M„) ,,J MH- (4-7) 

n=0 2 ^^2 + 

Here we explicitly write the projection \n){n\ for each species doubler. If one chooses the 
mass parameters so that 

Mo = -— < 0, M„ > for n ^ (4.8) 
a 

namely 

< mo < 2r (4.9) 
and if one lets all the mass parameters |M„| become large, one obtains 

"^^^ " ^^2[^ + 2M|] ^ = °' 

a75/^ (_l)"^,i[2 + i£-l^] /or n ^ 0. (4.10) 
If one chooses mo to satisfy 

2a|Mo| = 2mo = 1 (4.11) 

one recovers the correctly normalized continuum Dirac operator for the physical species 
and 75D ~ (— l)"75i for unphysical species doublers. In particular, the first relation in 
(4.10) can then be written as 

// = a75L> ~ 75ai ^ + 75(750^ ^)^ (4.12) 

which ensures the conventional Ginsparg- Wilson relation in the leading order. These 
properties become important in the following discussion. 

4.1 Generalized algebra with k = 1 

We now come back to the generalized algebra (2.1) with k = 1 

H^5 + ^5H^2H^ (4.13) 
10 



where H = a'j^D and = 75 — H^. This algebraic relation implies that 

'j,H' = [^,H + H^^]H - H[-f,H + H-f,] + H^^, = H^'y, (4.14) 
Namely, the algebraic relation (4.13) is equivalent to the two relations 

75i/2 - //%5 = 0. 
If one defines i^(3) = H^, the first relation of (4.15) becomes 

ff(3)75 + 75^^(3) =2fff3) (4.16) 

with Fs = 75 — i7(3), which is identical to the conventional Ginsparg- Wilson relation (1.1). 
We utilize this property to construct a solution to (4.15). Note that the operator F5 is 
identical in these three ways of writing in (4.13), (4.15), and (4.16). 

The physical condition for the operator H in (4.13) in the near continuum configuration 
is (Cf.(4.12)) 

~75ai^ + 75(75ai^)' (4.17) 
and thus if(3) in (4.16) should satisfy 

if(3) ~ [75az^ + 75(75a«^)^]3 

~ (75az^)3 + 75(75a^^)^ (4.18) 

as can be confirmed by noting 75 P+ ^75 = 0. Here only the leading terms in chiral 
symmetric and chiral symmetry breaking terms respectively are written. 
One can thus construct a solution for if(3) by 

^(3) = ^75 [1 + D'i) I ] (4.19) 



where we defined D^) byQ 



D^w^ ^ + (Bf - (4.20) 



The operators (^,B and the parameter mo/ a are the same as in the original Wilson fermion 
operator (4.1). By rewriting (4.19) as 

47.11+ 7.i/S'-^^] (4.21) 

in terms of the hermitian = 'J5D^} = [H^^Y and comparing it with (4.3), one 
can confirm that our operator if(3) satisfies the relation (4.16). The condition (4.18) is 
satisfied by noting 

Dp ^ + (Mf )3 (4.22) 

^It is also possible to use = + [B — ■^)'^: or any suitable (ultra-local) operator which 

satisfies 75i^[^^ = (75^1^V and (4.22). 
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in the near continuum configuration, where the mass parameters are given by 

(Mi3))3 , {ft3_(^)^(-)'-(^)^ft'-(^)^(^)'-(^« 

a a a a a a a a 

for ny^O. (4.23) 

Although we have the same condition on the parameters as before 

< mo < 2r (4.24) 

to avoid the species doublers, the value of mo itself is now required to satisfy 

2(mo)^ = 1 (4.25) 

to ensure the properly normalized physical condition (4.18). 

4.2 Reconstruction of H from i7(3) 

We now discuss how to reconstruct H, which satisfies (4.13), from H(2) defined above. 
The basic idea is to take a real cubic root of i7(3) as 

H^iH^-i))'/^ (4.26) 

in such a manner that H thus obtained satisfies the second constraint in (4.15). For 
this purpose, we first recall the essence of the general representation of the algebra (2.1) 
analyzed in Section 2, which is apphcable to (4.16) as well. 

If one defines the eigenvalue problem 

^(3)0n = {a\nf(f>n, (0n, 0n) = 1 (4.27) 

one can classify the eigenstates into the 3 classes: 
(i) n± ("zero modes"), 

^(3)0n = O, 75(/>„ = ±(/>n, (4.28) 

(n) N± ("highest states"), 

if(3)0„ = ±0„, 750„ = ±0„, respectively, (4.29) 
(iii) "paired states" with < |(aA„)^| < 1, 

H^s)<j)^ = {aXnf(pn, i/(3)(r50n) = -(aA,)3(r50n). (4.30) 

where 

r5 = 75-^(3). (4.31) 
12 



Note that r5(r50„) oc 0„ for < |(aA„)^| < 1. 
We obtain the index relation 

n 

= Xl(0n,750n) 
A„=0 

= n+ — n_ = index (4.32) 

where n± stand for the number of normahzable zero modes in the classification (i) above. 
We also have a chirality sum rule 

n+ + N+^n- + (4.33) 

where N± stand for the number of "highest states" in the classification (ii) above. 

If one denotes the number of states in the classification (iii) above by 2No, the total 
number of states (the dimension of the representation) N is given by 

N ^ 2{n+ + N+ + No) (4.34) 

which is expected to be common to all the fermion operators defined on the same lattice. 

Also, all the states 0„ with < |(aA„)^| < 1, which appear pairwise with (aXn)'^ ~ 
±|(aA„)^|, can be normalized to satisfy the relations 

l^(l>n = (aA„) + [1 - (aA„) Y^V-n, (4.35) 
where stands for the eigenstate with an eigenvalue opposite to that of 



Based on these general results in Section 2, we first observe that the index n+ — n_ in 
(4.32) is identical to the index of the expected solution of (4.13), although H(2,) satisfies 
(4.18). This observation is based on the relation 

n+ - n_ = ^(0„, T,f{{H^^)f/{aMf)4>n) (4.36) 

n 

which is valid for any regulator with /(O) — 1. One can perform the same analysis as in 
(3.7) in Section 3: The basic ingredient is the condition (4.18) for a physical momentum 
region in the smooth continuum limit and the absence of species doublers. The calculation 
analogous to (3.14) then gives 

n+ — n_ = lim Tr'^^fi—r) = lim Tr'y^gi-—r) (4.37) 

with g{x) = f{x^) and cjf(O) = 1. The right-hand side of this relation shows that the 
present index is identical to the index of the general operator in (2.1), which includes an 
expected solution of (4.13). Due to the chirality sum rule (4.33), we also obtain the same 
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value of — as for an expected solution of (4.13). 



The agreement of the index of if(3) with the index of the expected solution H of (4.13) 
suggests that we can define H operationally by 

H(t)n = aXn(j)n (4.38) 

by using the same set of eigenf unctions and (the cubic roots of) eigenvalues 

{04, {aXn} (4.39) 

as for if(3) in (4.27). Note that the operator Fs = 75 — i7(3) = 75 — H^, which reverses the 
signature of eigenvalues of "paired states" and defines the index, is consistently chosen to 
be identical for (4.16) and for (4.38)0. 

We can then confirm the second constraint in (4.15) and the defining algebraic relation 
(4.13) for any "paired state" (pn, 

[if%5 - 75i^^]0n = H^l5(pn " 75(aA„)Vn 

= H^{iaXnf(l)rr + [l-iaXnf]'^''(l).n} 

-(aA„)2{(aA„)30„ + [1 _ (aA„)T/'0-n} 

= (4.40) 

and 

[T,H + /fr5]0„ = r5(aA„)0„ - aA„(r50„) = (4.41) 

where we used the relations in (4.35) and the definition (4.38). For "zero modes" and 
the "highest states", which are the eigenstates of 75, the condition [-^^75 — 75-ff^]0„ = 
obviously holds, and the relation [T5H + HT^](f)n = is also confirmed. 

The general representation of the algebra (4.13) is obtained from the standard rep- 
resentation, which is defined by H in (4.38), 75 in (4.35), and the state vectors {4>n} in 
(4.39), by applying a suitable unitary transformation. 

5 Discussion 

When one considers the algebraic relation with a constant R 

75(75^) + (75/^)75 = 2Ra''^-'\j,D)"^^' (5.1) 
instead of (2.1), one can eliminate the paprameter i? by a scale transformation 

D^D' = (52) 

The path integral 

J VipVtpexp[J ^pD'^] (5.3) 



^This means that an explicit calculation of the chiral Jacobian (and chiral anomaly) for the theory 
defined by (4.13) is performed by TrT^ — Tr{'-f5 — i?(3)) in terms of iJ(3) in (4.19). 
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is equivalent to 

J V^pVipexplJ ipDip] (5.4) 

after absorbing the parameter into ijj, at least in a well regularized lattice path 

integral. Consequently, the parameter R and also the factor a'^'^'^^ do not have an intrinsic 
physical significance^ 

In contrast, the power of (75!?)^'^"'"^ in the right-hand side of (5.1) has an intrinsic 
physical meaning. One may recall the near continuum expressions (4.12) and (4.17) 

H ~ 7502^ + 75(7502^)^ for k = 0, 

H ~ 7502^ + 75(7502^)^ for k = l (5.5) 

respectively. The first terms in these expressions stand for the leading terms in chiral 
symmetric terms, and the second terms in these expressions stand for the leading terms 
in chiral symmetry breaking terms. This shows that one can improve the chiral symmetry^ 
by choosing a large parameter k. 

The Dirac operator for such a general value of k is constructed by rewriting (2.1) as a 
set of relations (see (4.14)) 

H'j, - j,H' = 0, (5.6) 

with H = a'j^D. The first of these relations (5.6) becomes identical to the ordinary 
Ginsparg- Wilson relation (1.1) if one defines H(^2k+i) = H"^^^^. One can construct a 
solution to (5.6) by following the prescription in Section 4 

where 

ct 

The operator H is then finally defined by (in the representation where H{2k+i) is diagonal) 

H = {H^^k+Df'^"^' (5.9) 

in such a manner that the second relation of (5.6) is satisfied. This condition is in deed 
satisfied as a generalization of (4.40) in the representation where H(2k+i) is diagonal. We 
use the relation (2.23) in this proof. Also the conditions < mg < 2r and 

2mf+^ = 1 (5.10) 

^ However, when one includes a Yukawa interaction, for example, this scaling argument need to be 
refined. 

^To avoid the mis-understanding, we note that the improvement of chiral symmetry here is meant 
in the sense of Wilsonian renormalization group. The chiral symmetry breaking term becomes more 
irrelevant for larger fc, and this should be interesting from a view point of regularization of field theory 
in general. Also, the approach to the continuum Dirac operator is controlled by two parameters, for 
example, by letting k — > large and a — > small simultaneously. 
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ensure a proper normalization of the Dirac operator H. 

However, one need to use a large enough lattice to accomodate the operator H with 
a large /c, since the operator (5.8) correlates lattice points far apart from each other for a 
large k. An explicit analysis of the locality property of our operator H as in Ref.[15] is left 
as an important problem. In the context of lattice simulation, it would be interesting to 
see how the chiral properties are modified if one uses the operator with k — 1, which has 
been analyzed in detail in this paper, instead of the conventional overlap Dirac operator 
with A; = 0. To detect the possible effects of A; 7^ in a reliable way, it is expected that 
one would have to consider a sufficiently large lattice and those observables which are 
sensitive to low energy excitations. 

As for the chiral fermions on the lattice, our general algebra (2.1) satisfies the decom- 
position 

(1 + 75)^(1^^ 

2 2 2 2 ^ ' 

with 

75 = 75 - 2a^^+\-i,Df^+\ {%f = 1 (5.12) 

by noting 75(75/^)^ = (751^)^75. This decomposition has the same form as for the overlap 
operator D satisfying the ordinary Ginsparg- Wilson relation. It is thus expected that one 
can apply the same considerations as in Refs.[16] and [17] to our general Dirac operator 
also. In particular, the fermion number non-conservation of the chiral theory defined by 

^ j Vi^Vi^eMj i^^^^D^^^i^} (5.13) 
follows from the fermion number transformation 

i) e*"'0, -0 ^ V'e"'". (5.14) 

If one remembers that the functional spaces of the variables ip and ip are specified by the 
projection operators (1 — 75)/2 and (1 + 75)/2, respectively, the Jacobian factor for the 
transformation (5.14) is given by [16] 

= exp{iaTr[75 - (750^)^*^+^]} = exp{m[n+ - n_]} (5.15) 

where the index is defined in (2.19). 

In conclusion, we have shown that the general idea of Ginsparg and Wilson can be 
precisely realized as a closed algebraic relation (2.1) and it admits of the explicit construc- 
tion of an infinite tower of new lattice Dirac operators as a generalization of the overlap 
Dirac operator. This should be interesting in the context of the regularization of field 
theory in general. 
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